Warm pasta phase in the Thomas-Fermi approximation 
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In the present paper, the pasta phase is studied at finite temperatures within a Thomas- Fermi (TF) 
approach. Relativistic mean field models, both with constant and density-dependent couplings, are 
used to describe this frustrated system. We compare the present results with previous ones obtained 
within a phase-coexistence description and conclude that the TF approximation gives rise to a richer 
inner pasta phase structure and the homogeneous matter appears at higher densities. Finally, the 
transition density calculated within TF is compared with the results for this quantity obtained with 
other methods. 
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I. INTRODUCTION 

A complete theoretical description of the processes in- 
volved in supernova explosion and stellar evolution de- 
pends on the construction of an adequate equation of 
state, able to describe matter ranging from very low den- 
sities to few times saturation density and from zero tem- 
perature to around 100 MeV. Protoneutron stars are be- 
lieved to be the initial phase of the remnants of a su- 
pernova explosion. They cool down to become a cold 
neutrino-free neutron star by emitting neutrinos, which 
means that the neutrino mean free path inside the star is 
an important quantity in the understanding of the stellar 
evolution. The neutrino mean free path depends on the 
reactions that may take place inside the star. Hence, the 
composition of the star plays a definite role in its fate. 

Protoneutron stars are believed to have a crust where 
a special matter known as pasta phase is expected to 
be present. The pasta phase is a frustrated system that 
arises in the competition between the strong and the elec- 
tromagnetic interactions [1-5] . The pasta phase appears 
at densities of the order of 0.001 - 0.1 fm~ 3 [5, 6] in 
neutral nuclear matter and in a smaller density range 
[7, 8] in /3-equilibrium stellar matter. The basic shapes 
of these structures were named [1] as droplets (bubbles) , 
rods (tubes) and slabs for three, two and one dimensions 
respectively. The ground-state configuration is the one 
that minimizes the free energy, i.e., the pasta phase is the 
ground-state configuration if its free energy per particle 
is lower than the one corresponding to the homogeneous 
phase at the same density. 

In previous works [6, 8] we have studied the existence 
of the pasta phase at zero and finite temperature within 
different parametrizations of the relativistic non-linear 
Walecka model (NLWM) and of the density-dependent 
hadronic model. In the present work we will consider 
the NLWM parametrization NL3 [9] and the density- 
dependent hadronic model TW [10]. These two models 
do not satisfy some of the constraints on the high den- 



sity behaviour of the equation of state (EOS) discussed in 
[11]. However, both models describe reasonably well the 
groundstate properties of stable and unstable nuclei, and, 
therefore, are adequate to study the EOS at subsatura- 
tion densities. Moreover, we are interested in studying 
the effect of the density dependence of the symmetry en- 
ergy on the pasta phase. Therefore, we have chosen a 
model (NL3) with a very large symmetry energy slope 
at saturation, above the limit imposed by isospin diffu- 
sion [12], and another one with a slope close to the lower 
limit defined by the same experiments. In both works 
[6, 8] two different methods were used: the coexisting 
phases (CP), both at zero and finite temperature, and the 
Thomas-Fermi (TF) approximation at zero temperature 
only. It was found that if /3-equilibrium is imposed the 
pasta phase does not appear in a CP calculation for the 
same surface energy parametrization used for fixed pro- 
ton fractions. This indicates the necessity to use a good 
parametrization for the surface energy which is temper- 
ature, proton fraction and geometry dependent, as also 
stressed in [13, 14]. The specific problem of an appropri- 
ate parametrization for the surface energy was tackled in 
[15]. 

In [16] it is found that the diffusion coefficients, related 
to neutrino opacities in dense matter, are significantly al- 
tered by the presence of the nuclear pasta in stellar mat- 
ter. These differences in neutrino opacities certainly in- 
fluence the Kelvin- Helmholtz phase of protoneutron stars 
as well as supernova explosion simulations. 

In the present paper we use the Thomas- Fermi approx- 
imation to obtain the pasta phase for two parametriza- 
tions, namely NL3 and TW, and compare the results 
with the ones obtained with the more naive coexisting- 
phases method at finite temperature. The CP method 
used here is improved with respect to references [6, 8] 
because it takes into account a better description of the 
surface energy, as in [15]. We check the differences in 
phase transition from pasta to homogeneous matter and 
see how different the internal structures of the pasta are. 
Some preliminary results for the pasta phase obtained at 
finite temperature within the Thomas Fermi approxima- 
tion have been published in [17, 18]. 

It was shown in [19] that the pressure and density at 
the inner boundary of the crust (transition pressure and 
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transition density) define the mass and moment of inertia 
of the crust. This establishes a relation between the equa- 
tion of state (EoS) and compact-stars observables. Sev- 
eral methods have been used to calculate this density for 
/3-equilibrium cold matter: a local equilibrium approxi- 
mation [20, 21], the thermodynamical spinodal method 
[6, 8, 21, 22], and the dynamical spinodal method used 
in [6, 8, 23] which predicts a transition density very close 
to the Thomas-Fermi result. It is expected that the tran- 
sition density lies in the metastable region between the 
binodal surface and the dynamical spinodal surface. We 
will compare the different predictions for the transition 
density at several temperatures and isospin asymmetries 
obtained from the thermodynamical spinodal, the dy- 
namical spinodal, the binodal surface and the TF and 
CP calculations. This will allow us to make an estima- 
tion of the error done within each method. 

In section II the basic expressions for the non-linear 
Walecka model are outlined. In section III A we briefly 
describe the TF approximation at finite temperature and 
in section IIIB the CP method is reproduced. In section 
IV the results are shown and discussed and in section V 
the final conclusions are drawn. 

II. THE FORMALISM 

We consider a system of protons and neutrons with 
mass M interacting with and through an isoscalar-scalar 
field <f> with mass m s , an isoscalar- vector field with 
mass m v , an isovector- vector field b M with mass m p . We 
also include a system of electrons with mass m e . The 
Lagrangian density reads: 

C= Ci + Ce+Ca+Cu+Cp+Cy, (1) 

i— p,n 

where the nucleon Lagrangian reads 

£i=&[7 M iD"-M*]Vi, (2) 

with 

iD» = id" - T V V» - ■ W - eli^A", (3) 
M* = M-T s (j). (4) 

The electron Lagrangian density is given by 

£ e = 4>e [7m + eA") - m e ] Ve, (5) 

and the meson Lagrangian densities are 

£ a = + l - (d^<f> - m 2 s <f> 2 - l -n4? - ^A0 4 ) 

r — —-F F^ v 

7 — 4 f*" ' 



where 0^ = d^-d^ , B M „ = d^h l/ -<9„b M -r p (b M x 
h v ) and F pv — d l _ L A l/ —d v A, x . The parameters of the mod- 
els are: the nucleon mass M — 939 MeV, the coupling 
parameters r s , r„, T p of the mesons to the nucleons, 
the electron ma ss m e a nd the electromagnetic coupling 
constant e = ^/47r/137. In the above Lagrangian den- 
sity r is the isospin operator. When the TW density 
dependent model [10] is used, the non-linear terms are 
not present and hence k = X = and the density depen- 
dent parameters are chosen as in [10, 24, 25]. When the 
NL3 parametrization is used, r, is replaced by gi, where 
i — s, v, p as in [9]. The bulk nuclear matter properties of 
the models we use in the present paper are given in Ta- 
ble I. We also include in the table some properties at the 
thermodynamical spinodal surface: p s is the upper bor- 
der density at the spinodal surface for symmetric matter 
(it defines the density for which the incompressibility is 
zero), p t and P t are, respectively, the density and the 
pressure at the crossing between the cold /3-equilibrium 
equation of state and the spinodal surface. They give a 
rough estimate of the density and pressure at the crust- 
core transition [6, 8]. The surface-tension coefficient o 
(see section III B) is also given for T = and two proton 
fractions, Y p — 0.5 and Y p — 0.3. 



TABLE I: Symmetric nuclear matter properties at the satura- 
tion density and at the spinodal surface. The surface tension 
a, the saturation density po and the corresponding binding 
energy B/A are also given for the proton fraction Y p = 0.3. 





NL3 


TW 




[9] 


[10] 


Po (fm- 3 ) 


0.148 


0.153 


B/A (MeV) 


16.3 


16.3 


po (Y p = 0.3) (fm" 3 ) 


0.118 


0.135 


B/A(Y P = 0.3) (MeV) 


10.9 


11.2 


K (MeV) 


272 


240 


£ sym . (MeV) 


37.4 


32.8 


M'/M 


0.60 


0.55 


L (MeV) 


118.3 


55.3 


K sym (MeV) 


100.5 


-124.7 


Qo (MeV) 


203 


-540 


K T (MeV) 


-698 


-332 


Ps (fm- 3 ) 


0.096 


0.096 


Pt (fm- 3 ) 


0.065 


0.085 


Pt (MeV/fm 3 ) 


0.396 


0.455 


a(Y p = 0.5) (MeV/fm 2 ) 


1.123 


1.217 


a(Y p = 0.3) (MeV/fm 2 ) 0.254 


0.613 



We next give the basic expressions for the construc- 
tion of the pasta phase within the Thomas-Fermi calcu- 
lation at finite temperature and of the coexisting-phases 
method. 
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III. THE PASTA PHASE 
A. Thomas-Fermi approximation 

In the present work we repeat the same numerical pre- 
scription given in [6] where, within the Thomas-Fermi 
approximation of the non-uniform npe matter, the fields 
are assumed to vary slowly so that the baryons can be 
treated as moving in locally constant fields at each point. 

From a formal point of view the Thomas-Fermi ap- 
proximation can be considered as the zeroth order term 
of a semiclassical expansion in the relativistic mean field 
theory derived within the framework of the relativistic 
Wigner transform of operators [26, 27]. Here, we take a 
more straightforward way to obtain the finite tempera- 
ture semiclassical TF approximation based on the density 
functional formalism. We begin with the grand canonical 
potential density: 

w = ^{/i+}, {/*-}, 0o,Vb,6o) = £t-TSt- ^2 ViPi , 

i—p.n,e 

(6) 

above i=p,n,e stands for the protons, neu- 

trons and electrons positive(negative) energy distribution 
functions and S t = S + S e , £ t = £ + £ e are the total en- 
ergy and entropy densities respectively. The total energy 
has been defined in [6] and for the entropy we take the 
one-body entropy: 

St = ~ £ /^/S {/i+(r,p)In/ j+ (r,p) (7) 

i—n,p,e 

+ [1 - / i+ (r, p)] In [1 - f i+ (r, p)] + (f i+ o /,_)} . 

where the ground-state (equilibrium) distribution func- 
tions are 



where 



/»±(r,p) 



/e±(r,p) 



1 



1 + exp [(e*(r,p)T^(r))/T] 



l + exp[(e e T^(r))/T]' 



,i=p,n 



(8) 
(9) 



with e*(r,p) = yV + M*(r) 2 , M*(r) = M — r s o (r), 

and e e = y/p 2 + m%. v e ( r ) = Pe + eA (r) is the elec- 
tron effective chemical potential and the nucleon effective 
chemical potential, i/j, i — p,n, is given by: 

v i = iH - T V V - y m b - |(1 + t 3i )A - (10) 

where the rearrangement term is given by [24]: 
^Rr ^ dT v , wn dT p 6 (r) dT s 



dp 



2 dp 



The equations of motion for the meson fields (see 
follow from the variational conditions: 

s s s 



(11) 



= d 3 ruj({f l+ }, {ft-}, 4> , V , b ) 



(12) 



The numerical algorithm for the description of the neu- 
tral npe matter at finite temperature is a generalization 
of the zero temperature case which was discussed in de- 
tail in [6] . The Poisson equation is always solved by using 
the appropriate Green function according to the spatial 
dimension of interest and the Klein-Gordon equations are 
solved by expanding the meson fields in a harmonic os- 
cillator basis with one, two or three dimensions based on 
the method proposed in [28] . The most important source 
of numerical problems are the Fermi integrals, hence, we 
have used an accurate and fast algorithm given in ref. [29] 
for their calculations. 



B. Coexisting phases 

In [6] a complete description of the coexisting-phases 
method applied to different parametrizations of the 
NLWM is given. In the description of the equations of 
state of a system, the required quantities are the bary- 
onic density p, the energy density £, the pressure P and 
the free energy density J 7 , explicitly written in [6, 8]. 

Two possibilities are discussed in the following: nuclear 
matter with fixed proton fraction and /?-equilibrium stel- 
lar matter. In the first case, electrons are related to the 
nuclear matter by the imposition of charge neutrality in 
such a way that the electron density is equal to the proton 
density. The equations 



v - El 



Pe = Pp 



(13) 



hold. In the second case, the fractions of nucleons and 
electrons are defined by the conditions of chemical equi- 
librium and charge neutrality. In this case, the enforced 
conditions are: 



Pp Pn Pe, Pe Pp 



(14) 



for neutrino-free matter and 



\r Pe^Pv /I c\ 

Pp = Pn -Pe + Pu, Yi = , p e = p p , (15) 



if trapped neutrinos are considered. 

As in [3, 6], for a given total density p, the pasta struc- 
tures are built with different geometrical forms in a back- 
ground nucleon gas. This is achieved by calculating from 
the Gibbs conditions the density and the proton fraction 
of the pasta and of the background gas. 

In building the pasta phase, the density of electrons is 
uniform. The total pressure is given by P = P 1 + P e + P v 
and the total energy density of the system is given by 



fo(r) Wo(v) Sbo(v)' 



£ = f£' + (1 - f)£ 11 +£ e +£ V + £surf + £ C o 



(16) 
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where I and II label the high and low density phase re- 
spectively and / is the volume fraction of phase I. Notice 
that matter with fixed proton fraction is neutrino-free 
and hence the neutrino pressure and energy density are 
zero. By minimizing the sum £ SU rf + £coui with respect 
to the size of the droplet /bubble, rod/tube or slab we get 
[3] £ aur f = 2£coui, and 

Seoul = jJ(eW /3 {°D{pl p»)) V * , (17) 

where F = f for droplets and F — 1 — / for bub- 
bles, <t is the surface energy coefficient, D is the di- 
mension of the system and the geometric factor $ = 

( 2 ~ D d- 2 2/D +F ) "DT2' D = X ' 2 ' 3 - Thc surface coeffi- 
cient a, necessary in the above expressions, is given by 
[30-32]: 

which is adequate for the parametrization of the depen- 
dence of a on the temperature. In the Appendix we show 
how to obtain this expression and the equivalence be- 
tween this expression and eq. (3.14) of reference [33]. 

According to the calculations performed in [15] the sur- 
face energy in terms of the proton fraction varies consid- 
erably from the NL3 to the TW parametrization. In thc 
first case it is always smaller. This quantity plays an im- 
portant role on the size and structures of the pasta phase 
discussed in the present paper. 

At this point it is worthy pointing out that the de- 
pendence of the energy density on the electromagnetic 
and surface contributions is commonly known as finite 
size effects. In [34-37] it was shown that for a weak sur- 
face tension the EoS obtained in a mixed phase resembles 
the one obtained with a Gibbs construction while, for a 
strong surface tension, thc Maxwell construction is re- 
produced. 

IV. RESULTS AND DISCUSSIONS 

In the present section we discuss the properties of the 
pasta phase obtained within thc two methods described 
above. We consider several temperatures, T — 0, 5, 7, 8 
MeV and two proton fractions Y p = 0.5 and 0.3. We 
only show results up to 8 MeV because it is not clear 
whether our framework model is realistic above this tem- 
perature , since thermal fluctuations are not taken into 
account. The problem of the effect of thermal fluctua- 
tions on the pasta structures has been studied by [38, 39] 
and it was shown that thermally induced displacements 
of the rod-like and slab-like nuclei can melt the lattice 
structure when these displacement are larger than the 
space available between thc cluster and the boundary of 
the Wigner-Seitz cell. While for the rod like clusters 
T = 8 MeV would still be acceptable, for the slabs T = 8 
MeV could already be too large. 



We perform most of our calculations for a fixed proton 
fraction, although matter in /3-equilibrium with trapped 
neutrinos is known to be important. In fact during 
the first seconds of the protoneutron star, neutrinos are 
trapped. They start to diffuse out 10-15 s after the super- 
nova explosion [40, 41], when the fraction of leptons de- 
creases from a maximum value, ~ 0.4, constant through 
out the star, to smaller values. Neutrino emission is es- 
sential to understand neutron star cooling [42, 43]. How- 
ever, neutrinos couple weakly to nuclear matter, therefore 
for the pasta calculation their presence mainly defines the 
proton fraction of the warm stellar matter. In [23], we 
have shown for a wide set of models, including NL3 and 
TW, that for the density range of interest for the pasta 
phase, the proton fraction for the largest trapped frac- 
tion of neutrinos, which corresponds to a lepton fraction 
Yi ~ 0.4, is approximately 0.3. Moreover, in a recent pa- 
per [15] we have studied the pasta phase in /3-equilibrium 
stellar matter with trapped neutrinos (V; = 0.4) within 
the more schematic coexisting phases approach and we 
have seen that the pasta phase structure and its exten- 
sion is very similar to the one obtained for a fixed proton 
fraction of about 0.3. 

Since the inclusion of trapped neutrinos in the calcu- 
lation does not bring much more information about the 
pasta phase itself (the main point of the present work), 
we have only included neutrinos in the present calculation 
for two cases (NL3 and TW, T = 5 MeV and V; = 0.4) 
in order to show the similarity with the Y p — 0.3 results. 

In Tables II and III we display the onset densities for 
each of the pasta structures and the homogeneous phase 
for the parametrizations NL3 and TW respectively. The 
pressures identifying the starting (P 1 ) and ending (P 2 ) 
points of the pasta phase are also given. We investigate 
all possible structures for T = 0,5,7 and 8 MeV and 
proton fraction equal to 0.5 (symmetric matter) and Y p = 
0.3. Both CP and TF approximations were used. 

From Table II one can observe that TF always presents 
a richer inner pasta phase structure as compared with the 
CP method. For T — 8 MeV and symmetric matter, the 
pasta phase no longer exists if the CP method is used, 
but almost all internal structures are present within the 
TF calculation. Generally, the homogeneous phase be- 
comes the ground-state matter for densities much higher 
within the TF approach than within thc CP method. 
Just some exceptions were found for a complete structure 
(3D, 2D, ID, 2D, 3D) with the TF calculation: T = 7 
MeV for Y p = 0.3 and T — 8 MeV for Y p = 0.3 and 
Y p = 0.5, which lack the slab configuration. In all other 
cases, all possible configurations were found. This is not 
true within a CP method, where many structures were 
missing, giving rise to a much simpler pasta phase. 

From Table III one sees that most of the conclusions 
drawn for the NL3 parametrization hold true also for 
TW, as a richer inner pasta phase structure obtained 
within the TF approach than within thc CP method. 
Once more the homogeneous phase becomes the ground- 
state matter for densities higher within the TF approach 
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FIG. 1: Pasta phases: comparison between NL3 (dashed line) and TW (full line) and the methods CP and TF for a) /3- 
equilibrium cold stellar matter; b) Y v — 0.5, c) Y v — 0.3. The thick lines refer to the Thomas-Fermi (TF) calculation and the 
thin red lines to the coexisting-phases (CP) calculation. 



than within the CP method. When the TW density 
dependent model is used, one can see that all possible 
internal pasta configurations are always found with the 
TF approximation. This is not the case with the NL3 
parametrization. Different pasta phase internal struc- 
tures lead to different diffusion coefficients causing dif- 
ferent neutrino mean free paths, as seen in [16]. These 
differences in neutrino opacities may have consequences 
in calculations of protoneutron star evolution. However, 
although the internal structures are different as compared 
with the NL3 parametrization, the densities at which the 
pasta structure appears and ends are quite similar. 

In tables II and III we have also included the results for 
/3-equilibrium matter with trapped neutrinos and a lep- 
ton fraction Yi = 0.4 at T = 5 MeV. It is seen that these 
results are very similar to the ones obtained for Y p = 0.3. 
This is understandable because, for the pasta phase den- 
sities, the proton fraction changes from ~ 0.29 for the 
lowest densities to ~ 0.32 for the largest ones. A proton 
fraction smaller than 0.3 at the onset of the pasta phase 
explains the onset of the pasta at slightly larger densi- 



ties. On the other hand a proton fraction Y p > 0.3 at 
the upper limit of the pasta phase corresponds in the TF 
calculation to a slightly larger transition density than the 
Y p = 0.3 value, and lies between the value obtained for 
Y p = 0.5 and Y p = 0.3. In the CP calculation the result 
also lies between Y v = 0.5 and Y p = 0.3 but is instead 
smaller than the Y p — 0.3 value because CP does not 
allow for the rearrangement of the proton distribution. 

The comparison between Tables II and III is more 
easily done analysing Fig. I. In Fig. 1 we compare 
the density range for which each pasta configuration 
exists within NL3 and TW. The thick lines stand for 
the Thomas-Fermi calculation, the thin ones for the CP 
method. Full lines represent TW and dahsed ones NL3. 

For symmetric matter the main difference is the the 
appearance of the different phases at slightly smaller 
densities within NL3. The largest differences occur for 
Y p = 0.3: NL3 has no slab phase at T=7 and 8 MeV and 
quite large tube and bubble phases. 

It has been discussed in [5, 22] that the characteristics 
of the pasta phase are strongly related with the density 
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dependence of the symmetry energy. In particular, in 
[5] the pasta phases have been calculated using quantum 
molecular dynamics. The authors have considered two 
models, and have obtained a quite different structure for 
the pasta phase of both models. The slab phase is very 
small in one of the models and although larger in the 
other, it is one of the phases that first disappears when 
the temperature is raised. 

A large value of the slope of the symmetry energy L 
means a smaller symmetry energy at sub-saturation den- 
sities. Therefore, since TW has a much smaller symmetry 
energy slope, L = 55 MeV for TW and L = 118 MeV for 
NL3, at sub-saturation densities its symmetry energy is 
larger and the neutron gas equation of state has a larger 
energy. As a consequence, neutrons do not drip so eas- 
ily, the surface tension is larger and the surface is less 
diffuse. In Table I we have included the value of the sur- 
face tension a of TW and NL3 for the proton fractions 
0.5 and 0.3. It is seen that TW has a much larger a 
for Y p = 0.3 than NL3 which is in accordance with its 
symmetry energy slope. We may now understand the 
different properties of the TW pasta: the larger surface 
energy explains the fact that all the pasta phases ex- 
tend to larger densities. For non-symmetric matter as 
for Y p = 0.3 this effect is even stronger. This explains 
some of the largest differences between NL3 and TW as 
pointed out before, namely the fact that NL3 has no slab 
phase above T—7 MeV. As already noticed in [6, 36] the 
extension of the pasta phase decreases with the increase 
of the temperature. This is expected because the surface 
tension decreases with temperature. 

Fig. 1 also allows for a clear comparison between the 
CP and TF methods: the phases droplets, rods and slabs 
are stable until larger densities with CP. On the other 
hand, the phases tubes and bubbles do not appear as the 
ground-state configuration since they have a higher free 
energy than the homogeneous matter. Other important 
differences are the smaller extension of the pasta phase 
within CP (it starts at larger densities and finishes at 
smaller densities), and the disappearance of the pasta 
phase at smaller temperatures. A smaller extension of 
the pasta phase within the CP method is partially due 
to non-selfconsistent treatment of the Coulomb interac- 
tion [3], which prevents the rearrangement of the pro- 
ton distributions (see also [15]). This is more strongly 
seen for the Y p — 0.5 when Debye screening effects are 
stronger. We have calculated the critical temperatures 
within the CP calculation, above which the pasta phase 
does not exist. For NL3 we have T C CP = 7.73, (8.09) 
MeV for Y p = 0.5, (0.3), for TW T C CP = 7.95, (8.42) 
MeV for Y p = 0.5, (0.3). The CP critical temperature 
occurs when the free energy of the homogeneous phase 
is smaller than the free energy of the pasta phase, and is 
not defined by a zero surface tension, as in TF. At and 
above the critical temperature, pasta clusters still exist 
within CP calculation but with a larger free energy than 
homogeneous matter. 

Finally, it is interesting to compare the prediction of 



different methods for the crust-core (non-homogeneous- 
homogeneous) transition density. The thermodynamical 
binodal and spinodal of proton-neutron (pn) matter give 
a good estimation of the extension of the pasta phase 
in nuclear matter. The binodal surface is defined in the 
p,Y p ,T space as the surface where the gas and liquid 
phases coexist and it defines an upper limit for the ex- 
tension of the pasta phase since it does not take into 
account neither Coulomb nor finite size effects. The ther- 
modynamical spinodal defines the surface in the p, Y p , T 
space where the curvature of the free-energy of pn mat- 
ter becomes negative, and gives a rough estimation of the 
lower limit of the pasta phase extension. Also in this case 
nor electron effects neither finite range effects are taken 
into account. The transition density is obtained from 
the intersection between the /3-equilibrium EoS and the 
spinodal or binodal surfaces. It is possible to get a better 
estimation of the lower limit of the pasta phase extension 
if, instead of the thermodynamical spinodal, the dynam- 
ical spinodal is calculated. This surface is defined by the 
density, proton fraction and temperature for which mat- 
ter becomes unstable to small density fluctuations and is 
calculated including electrons, the Coulomb interaction 
and finite range effects. It has been shown in [6, 8, 21, 22] 
that the thermodynamical method gives a good estima- 
tion of the transition density for cold /3-equilibrium mat- 
ter, although a bit too large. For a fixed temperature, 
the binodal, thermodynamical and dynamical spinodals 
touch each other (or almost touch) at the critical point 
[44]: this corresponds to a quite small proton fraction 
and is close to the crust-core transition density. How- 
ever, for less asymmetric matter as the one that occurs 
in stellar matter with trapped neutrinos, the binodal oc- 
curs at quite larger densities than the spinodals. This is 
clearly seen in Fig. 2 comparing the beta-equilibrium re- 
sults for cold stellar matter with the results for Y p = 0.5 
or 0.3. 

For densities above the binodal surface, matter is ho- 
mogeneous. For densities below the dynamical spinodal 
matter is non-homogeneous. Between this two surfaces 
we may find matter in a metastable configuration. The 
most probable configuration is the one with the smallest 
free energy density. The results given in Tables II and III 
refer to the energetically favoured configurations. How- 
ever, we may argue that at finite temperature there may 
be a coexistence of several configurations with larger or 
smaller probability according to the corresponding free 
energy density. This partially justifies the large differ- 
ences between the TF and CP results: although the 
ground-state configurations are different, there is not a 
large energy difference between the configurations and, 
taking into account an average over all the possible con- 
figurations with a correct probability factor would give 
similar predictions. 

In Fig. 2 we compare the transition density from a non- 
homogeneous phase to a homogeneous phase obtained 
from the binodal surface [44], the dynamical spinodal 
surface [20, 21, 23, 45-47], the thermodynamic spinodal 
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FIG. 2: Transition density, for several temperatures T = 0, 5, 7 and 8 MeV and proton fractions Y p = 0.5 (upper plots), 
0.3 (lower plots) and /3-equilibrium (eq-/3) at T=0 (left part of upper plots), obtained using different methods within the 
models a) NL3 and b) TW. The methods used are: estimation from binodal (Bin), estimation from dynamical spinodal (Sp-d) 
estimation from thermodynamical spinodal (Sp-th), calculation using the coexisting- phases method (CP) and calculation using 
a Thomas- Fermi method (TF). 



surface [21, 48, 49] and from the two methods we have 
used for the calculation of the pasta phase. As expected 
the TF result lies always between the results obtained 
from the dynamical spinodal and the binodal surfaces. 
This is a self-consistent method that should satisfy these 
two constraints. This is no longer true for the CP cal- 
culation, which, for larger temperatures, fails with pre- 
dictions that lie below the dynamical spinodal result. It 
is impressive that the thermodynamical-spinodal result 
gives so good results when it does not take into account 
neither the surface nor the Coulomb effects. These are 
good news because it allows a quite safe prediction from 
a simple calculation. Finally, it should also be pointed 
out that for /3-equilibrium matter all methods except CP 
give similar results. This is due to the occurrence of the 
transition density close the critical point where the spin- 
odal and binodal surfaces touch and the pressure on these 
surfaces is maximum. 



V. CONCLUSIONS 

In the present work we have calculated the pasta phase 
at zero and finite temperature applying two different 
methods already used in previous works: the more naive 
coexisting-phases method based on the Gibbs construc- 
tion and the self-consistent Thomas Fermi method. For 
the first method we need to know the surface tension of 
the models as a function of temperature and proton frac- 
tion. This quantity was parametrized from a Thomas 
Fermi calculation for semi-infinite nuclear matter [15]. 
We had already compared the two models at zero tem- 



perature [6, 8] but for the CP calculation we had used 
a parametrization for the surface tension obtained from 
Skyrme forces. Since the appearance of the pasta phases 
has a strong dependence on this quantity the comparison 
between the models at zero temperature was not conclu- 
sive. 

We have considered two different relativistic nuclear 
models: NL3, a parametrization of the NLWM, and TW, 
a parametrization of the density dependent hadronic 
model. These two models have quite different behaviors 
both in the isoscalar and isovector channels. In particu- 
lar, TW has a quite small slope of the symmetry energy 
at saturation, L — 55 MeV, while this value is 118 MeV 
for NL3. The properties of the models are clearly re- 
flected in the pasta phase structure as discussed in [5] . 

We conclude that while the CP method allows the 
determination of the overall trends of both models, it 
fails at a more detailed and quantitative level. The 
main trends observed are: the model having a larger sur- 
face tension predicts larger density ranges for the pasta 
phases; the pasta phase extension decreases with tem- 
perature; some of the less stable pasta structures may 
disappear for larger temperatures and/or isospin asym- 
metries. We have also shown that the overall conclusions 
obtained within the Thomas Fermi calculation at finite 
temperature agree with the conclusions obtained from a 
quantum molecular dynamics calculation [5]. In fact, it 
was shown that the structure of the pasta is sensitive to 
the model, namely the density dependence of the sym- 
metry energy. A model with a large symmetry energy 
slope such as NL3, has a quite small symmetry energy at 
sub-saturation densities and this favours neutron drip. 
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TABLE II: Pasta phases within NL3 for several temperatures and isospin asymmetries. The onset densities of the various 
phases (droplets, rods, slabs, tubes and bubbles) are shown together with the pressure of the homogeneous phase at the lower 
(P 1 ) and upper (P 2 ) border of the pasta phase. The results have been obtained within the coexisting phases method (CP) and 
the Thomas Fermi approximation (TF). 











onset density (fm 3 ) 




pi 1 p2 






EOS droplet 


rod slab tube 


bubble 


horn 


(MeV fm" 3 ) 


T 


= MeV 














Y P 


= 0.5 


















TF 


0.000 


0.025 0.043 0.073 


0.089 


0.114 


- / 2.75 






CP 


0.000 


0.032 0.053 - 




0.093 


- / 1.88 


Y p 


= 0.3 


















TF 


0.000 


0.019 0.036 0.059 




0.099 


- / 1-11 






CP 


0.000 


0.030 0.048 0.084 


- 


0.093 


- / 0.96 


P- 


equilibrium 


TF 
CP 


0.000 


- 

- - - 


- 
- 


0.055 


- /0.26 


T 


= 5 MeV 














y P 


= 0.5 


















TF 


0.0020 


0.026 0.043 0.067 


0.081 


0.106 


0.024 / 2.55 






CP 


0.0035 


0.031 0.051 - 


- 


0.078 


0.046 / 1.53 




= 0.3 


















TF 


0.0022 


0.024 0.038 0.053 


0.066 


0.092 


0.018 / 1.08 






CP 


0.0029 


0.032 0.049 0.081 


- 


0.082 


0.025 / 0.85 


Yi 


= 0.4 


















TF 


0.0024 


0.023 0.036 0.054 


0.065 


0.093 


0.023 / 5.24 






CP 


0.0033 


0.031 0.049 - 


- 


0.080 


0.024 / 1.43 


T 


= 7 MeV 














Y P 


= 0.5 


















TF 


0.0048 


0.026 0.042 0.054 


0.072 


0.098 


0.076 / 2.37 






CP 


0.013 


0.030 0.049 - 


- 


0.056 


0.22 / 1.06 


Y P 


= 0.3 


















TF 


0.0052 


0.025 - 0.036 


0.059 


0.084 


0.056 / 1.02 






CP 


0.0107 


0.032 0.048 - 




0.065 


0.112 / 0.67 


T 


= 8 MeV 














Y P 


= 0.5 


















TF 


0.0069 


0.025 - 0.036 


0.063 


0.092 


0.121 / 2.23 






CP 










-/- 


Y P 


= 0.3 


















TF 


0.0076 


0.026 - 0.033 


0.055 


0.078 


0.090 / 0.96 






CP 


0.0239 


0.032 - 




0.043 


0.253 / 0.44 



As a consequence the background neutron gas is larger 
for NL3 at a given density and therefore, the favoured 
pasta structures change shape at smaller densities. 

We have also compared the crust-core transition den- 
sity, from a non-homogeneous phase to a homogeneous 
phase, obtained from the above two methods as well as 
using other methods often referred to in the literature. 
It was also shown that the dynamical spinodal defines 
a lower limit while the binodal a higher limit and the 
TF result lies between the two limits. The CP methods 



fails these constraints both for large temperatures, above 
T=5 MeV, and very asymmetric matter. We have ob- 
tained the interesting result that the estimates obtained 
from a thermodynamical calculation are very close to the 
prediction of the TF calculation. 



9 



ACKNOWLEDGMENTS by FCT and FEDER (Portugal) under the projects 

PTDC/FIS/64707/2006, CERN/FP/109316/2009 and 
This work was partially supported by Capes/FCT SFRH/BPD/64405/2009 and by Compstar, an ESF Re- 
n. 232/09 bilateral collaboration, by CNPq (Brazil), search Networking Programme. 



D. G. Ravenhall, C.J. Pethick and J.R. Wilson, Phys. 
Rev. Lett. 50, 2066 (1983). 

C.J. Horowitz, M.A. Perez-Garcia and J. Piekarewicz, 
Phys. Rev. C 69, 045804 (2004); C.J. Horowitz, M.A. 
Perez-Garcia, D.K. Berry and J. Piekarewicz, Phys. Rev. 
C 72, 035801 (2005). 

T. Maruyama, T. Tatsumi, D.N. Voskresensky, T. Tani- 
gawa and S. Chiba, Phys. Rev. C 72, 015802 (2005). 
Gentaro Watanabe, Toshiki Maruyama, Katsuhiko Sato, 
Kenji Yasuoka, and Toshikazu Ebisuzaki, Phys. Rev. 
Lett. 94, 031101 (2005). 

Hidetaka Sonoda, Gentaro Watanabe, Katsuhiko Sato, 
Kenji Yasuoka, and Toshikazu Ebisuzaki, Phys. Rev. C 
77, 035806 (2008); ibid, Phys. Rev. C 81, 049902 (2010). 
S.S. Avancini, D.P. Menezes, M.D. Alloy, J.R. Marinelli, 
M.M.W. Moraes and C. Providencia, Phys. Rev. C 78, 
015802 (2008). 

J. Xu, L.W. Chen, B.A. Li and H.R. Ma, Phys. Rev. C 
79, 035802 (2009). 

S.S. Avancini, L. Brito, J.R. Marinelli, D.P. Menezes, 
M.M.W. de Moraes, C. Providencia and A.M. Santos - 
Phys. Rev. C 79, 035804 (2009). 

G. A. Lalazissis, J. Konig and P. Ring, Phys. Rev. C 55, 
540 (1997). 

S. Typel and H. H. Wolter, Nucl. Phys. A656, 331 
(1999); Guo Hua, Liu Bo and M. Di Toro, Phys. Rev. 
C 62, 035203(2000). 

T. Klahn et al, Phys. Rev. C 74, 035802 (2006). 
L.W.Chen, C. M. Ko, and B. A. Li, Phys. Rev. Lett. 94, 
032701 (2005); Phys. Rev. C 72, 064309 (2005). 
P. Gogelein, E.N.E. van Dalen, C. Fuchs and H. Miither, 
Phys. Rev. C 77, 025802 (2008). 

E. N.E. van Dalen, C. Fuchs and A. Faessler, Eur. Phys. 
J.A. 31, 29(2007). 

S.S. Avancini, C.C. Barros, D.P. Menezes and C. 
Providencia, Phys. Rev. C 82, 025808 (2010). 
M.D. Alloy and D.P. Menezes, arXiv:1011.0968 [nucl-th]. 
S. S Avancini et al. Int. J. Mod. Phys. D 19, 1587 (2010). 
Toshiki Maruyama, Toshitaka Tatsumi, 

arXiv:1009.1198vl [nucl-th]. 

Bennett Link, Richard I. Epstein, and James M. Lat- 
timer, Phys. Rev. Lett. 83, 3362 (1999). 
C. J. Pethick, D. G. Ravenhall, and C. P. Lorenz, Nucl. 
Phys. A 584, 675 (1995). 

B. A. Li, L. W. Chen, C. M. Ko, Phys. Rep. 464, 113 
(2008); J. Xu et al, Astrophys. J. 697, 1549 (2009). 
Kazuhiro Oyamatsu and Kei Iida, Phys. Rev. C 75, 
015801 (2007). 

C. Ducoin, C. Providencia, A. M. Santos, L. Brito, and 
Ph. Chomaz, Phys. Rev. C 78, 055801 (2008). 

T. Gaitanos, M. Di Toro, S. Typel, V. Baran, C. Fuchs, 
V. Greco and H. H. Wolter, Nucl. Phys. A 732, 24 (2004). 
S.S. Avancini, L. Brito, D. P. Menezes and C. 
Providencia, Phys. Rev. C 70, 015203 (2004). 
P. Ring and P. Schuck, The nuclear many-body problem, 



Springer- Verlag, 1980. 
[27] M. Centelles, X. Vinas and M. Barranco, Ann. Phys. 221 
165 (1993). 

[28] Y.K. Gambhir, P. Ring and A. Thimet, Ann. Phys. 198, 
132 (1990). 

[29] J. M. Aparicio, Astrophys. J. Suppl. Ser. 117, 627 (1998). 
[30] M. Nielsen and J. da Providencia, J. Phys. G: Nucl. Part. 

Phys. 16, 649 (1990). 
[31] C. Providencia, L. Brito and J. da Providencia, M. 

Nielsen and X. Vinas, Phys. Rev. C 54, 2525 (1996). 
[32] D. P. Menezes and C. Providencia, Phys. Rev. C 60, 

024313 (1999). 

[33] M. Centelles, M. Del Estal and X. Vinas, Nucl. Phys. A 

635, 193 (1998). 
[34] D.N. Voskresensky, M. Yasuhira and T. Tatsumi, Phys. 

Lett. B 541, 93 (2002); Nucl. Phys. A 723 291 (2003). 
[35] N. Yasutake and K. Kashiwa, Phys. Rev. D 79, 043012 

(2009). 

[36] T. Maruyama, N. Yasutake and T. Tatsumi, 

arXiv:1009.1196vl [nucl-th]. 
[37] T. Maruyama, S. Chiba, H. J. Schulze, T. Tatsumi, Phys. 

Rev. D 76 (2007) 123015. 
[38] C. J. Pethick and A. Y. Potekhin, Phys. Lett. B 427, 7 

(1998). 

[39] G. Watanabe, K. Iida and K. Sato, Nucl. Phys. A 676, 
455 (2000); Nucl. Phys. A 726, 357 (2003) . 

[40] A. Burrows and J. M. Lattimer, Astrophys. J. 307, 178 
(1986). 

[41] M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis, J. M. 

Lattimer, and R. Knorren, Phys. Rep. 280, 1 (1997). 
[42] G. D. Yakovlev, A. D. Kaminker, O. Y. Gnedin, and P. 

Haensel, Phys. Rep. 354, 1 (2001). 
[43] A. B. Migdal, EE. Saperstein, M.A. Troitsky, D.N. 

Voskresensky, Phys. Rep. 192, 179 (1990). 
[44] H. Miiller and B. D. Serot, Phys. Rev. C 52, 2072 (1995). 
[45] F. Douchin and P. Haensel, Phys. Lett. B 485, 107 

(2000). 

[46] C. Providencia, L. Brito, S. S. Avancini, D. P. Menezes, 
and Ph. Chomaz, Phys. Rev. C 73, 025805 (2006); L. 
Brito, C. Providencia, A. M. Santos, S. S. Avancini, D. 
P. Menezes, and Ph. Chomaz, Phys. Rev. C 74, 045801 
(2006); H. Pais, A. Santos, L. Brito, and C. Providencia, 
Phys. Rev. C 82, 025801 (2010). 

[47] C. Ducoin, J. Margueron, and P. Chomaz, Nucl. Phys. A 
809, 30 (2008). 

[48] B.-A. Li, A. T. Sustich, M. Tilley, and B. Zhang, Nucl. 

Phys. A 699, 493 (2002). 
[49] S. S. Avancini, L. Brito, D. P. Menezes, and C. 

Providencia Phys. Rev. C 70, 015203 (2004); S. S. 

Avancini, L. Brito, Ph. Chomaz, D. P. Menezes, and C. 

Providencia, Phys. Rev. C 74, 024317 (2006). 



10 



Appendix 

In this appendix we calculate Eq. (18) for the surface 
tension a, following very closely reference [32]. The sys- 
tem is composed by np matter and the density depends 
only upon the z coordinate. Notice that no Coulomb 
field is included in the calculation of a. We start from 
the grand-canonical potential, Eq. (12), where 



ft : 

with 



JdV jl [(V0 O ) 2 - (Wo) 2 - (V6 ) 2 ] - Ve/j 

(19) 



Vef = 



™ 2 ^o + + | A ^0 



+ 2T £ J (2.)3 

t—p,n 

+ ln(l + e-( £ *+^/ T ) 



«OT " m$ + 2E*p] 

hn(l + e-^'-^/ T ) 



(20) 



The equations of motion for the the meson fields are ob- 
tained by minimizing ft with respect to each field, 



V 2 ^o 
V 2 V 
V 2 fe 



dV, 



ef 



dz 2 
d 2 V 



dV f 



ef 



dz 2 dV 
d 2 b _ dV ef 



dz 2 



Using the relation 

d 2 W _ d_ ( dW 
dz 2 dz \ dz 



db 



1 d fdW 



(21) 
(22) 
(23) 



<j>o,V 


,b , 


we o 


d 


'\( 


' d<po 


dcj)o 




v dZ 


d 


~\( 


'dV 


Wo 




, dz 


d 




Ubo 


dbo 


\ 


\ dz 



2 dW V dz 
n 

-Vef 
'Vef 
-Vf 



(24) 



Summing adequately the three equations (the second and 
third equations multiplied by -1), we get 



= 0. (25) 



This is equivalent to saying that 

2 /,, r \2 /,, \ 2' 



d<j) 
dz 



dVp 

dz 



db 

dz 



+ V ef =C, (26) 



where C is a constant which corresponds to the bulk 
contribution to the grand canonical potential density and 
can be identified with the pressure P. Replacing Eq. (26) 
into (19) we obtain 



ft = / dV 



d(f>o 
dz 



dVp 

dz 



db Q 

dz 



CV. 



(27) 

The surface energy is obtained from the free energy of a 
system with a fixed number of particles N = N p + N n , in 
which a cluster of arbitrary size exists in the background 
of the vapor phase. The free energy reads 



= ft + fj,pN p + p n N n , 

S l-oo dZ ( dz ) 
-CV + n P N p + p n N n . 



dz 



dbo 

dz 



For a cluster of volume V and surface S, we have 

F = Sa - CV + HpNp + p n N n . 
The surface energy per unit area of this cluster is 



I-co ^ ( 



dV \ 2 fdb xJ 



dz 



dz 



(28) 



(29) 



(30) 



Eq. (29) can be rewritten in the form of Eq. (3.14) of 
reference [33] 



/oo 
dz [T{z) -T g - fJ,p(p p (z) - p p ,g) 
-OO 

-p n {p n {z) - p n ,g)] 



(31) 



where T g = -P + p p p p . g + p n pn, g , P P , g and p n<g are 
the free energy density, the proton density and the neu- 
tron density of the gas. We have checked numerically the 
equivalence between eqs. (30) and (31). 
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TABLE III: Pasta phases within TW for several temperatures and isospin asymmetries. The onset densities of the various 
phases (droplets, rods, slabs, tubes and bubbles) are shown together with the pressure of the homogeneous phase at the lower 
(P 1 ) and upper (P 2 ) border of the pasta phase. The results have been obtained within the coexisting phases method (CP) and 
the Thomas Fermi approximation (TF). 











onset density (fm 3 ) 




pi / p2 






EOS droplet 


rod slab tube 


bubble 


horn 


(MeV far 3 ) 


T 


= MeV 














y P 


= 0.5 


















TF 


0.000 


0.025 0.043 0.075 


0.090 


0.115 


- /2.77 






CP 


0.000 


0.033 0.054 - 


_ 


0.093 


- /1.85 




= 0.3 


















TF 


0.000 


0.021 0.040 0.071 


0.085 


0.109 


-/ 1.22 






CP 


0.000 


0.031 0.051 0.89 


_ 


0.099 


-/ 0.96 


/3-equilibrium 


















TF 


0.000 


_ _ _ 


_ 


0.073 


- / 0.36 






CP 


0.003 


_ _ _ 


_ 


0.066 


0.007/ 0.29 


T 


= 5 MeV 














Y P 


= 0.5 


















TF 


0.0022 


0.026 0.043 0.071 


0.084 


0.107 


0.027 / 2.53 






CP 


0.0038 


0.032 0.052 - 


_ 


0.079 


0.050 / 1.55 


Y P 


= 0.3 


















TF 


0.0024 


0.025 0.042 0.068 


0.080 


0.103 


0.020 / 1.20 






CP 


0.0032 


0.033 0.052 - 


- 


0.087 


0.028 / 0.86 


Yi 


= 0.4 


















TF 


0.0026 


0.025 0.042 0.068 


0.081 


0.104 


0.026 / 1.64 






CP 


0.0035 


0.033 0.051 - 




0.085 


0.036 / 1.13 


T 


= 7 MeV 














Y P 


= 0.5 


















TF 


0.0051 


0.027 0.043 0.063 


0.077 


0.100 


0.082 / 2.41 






CP 


0.0130 


0.031 0.051 - 




0.060 


0.225 / 1.16 


Y P 


= 0.3 


















TF 


0.0055 


0.027 0.042 0.062 


0.074 


0.096 


0.061 / 1.16 






CP 


0.0110 


0.033 0.051 - 




0.071 


0.120 / 0.72 


T 


= 8 MeV 














Y P 


= 0.5 


















TF 


0.0073 


0.027 0.043 0.056 


0.072 


0.095 


0.130 / 2.32 






CP 










-/- 


Y P 


= 0.3 


















TF 


0.0079 


0.030 0.043 0.057 


0.069 


0.091 


0.096 / 1.13 






CP 


0.0207 


0.033 0.050 - 




0.055 


0.234 / 0.57 
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Erratum: Warm pasta phase in the Thomas-Fermi 
approximation [Phys. Rev. C 82, 055807 (2010)] 



PACS number(s): 21.65.-f, 24.10.Jv, 26.60.-c, 95.30.Tg 



There was an error in the numerical code for the CP 
method that resulted in using the zero temperature sur- 
face tension for all the cases, instead of the temperature 
dependent one. Also, applying the CP method, we found 
that in the surface energy parametrization given in Ref. 
[1] the global proton fraction should be used. As a result, 
the agreement between the more "naive" CP and the TF 



calculations has improved. The figures and the table in- 
cluding the correct values for CP are shown below. 

In the following, some specific comments we made are 
corrected. In the CP approximation, the "pasta" phase 
is found in all the cases considered: the critical tempera- 
ture for its existence is therefore higher than 8 MeV. The 
region in which the "pasta" phase exists is still generally 
smaller in the CP approach than in the TF one, but some 
exceptions are found (see Fig. 2). Concerning the tran- 
sition density (see Fig. 1), the CP results are now much 
closer to the TF ones. Moreover, for temperatures lower 
than 7 MeV, they satisfy the constraints imposed by the 
two spinodals. 



[1] S.S. Avancini, C.C. Barros, D.P. Menezes and C. 
Providencia, Phys. Rev. C 82, 025808 (2010). 
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FIG. 1: Pasta phases: comparison between NL3 (dashed line) and TW (full line) and the methods CP and TF for a) /3- 
equilibrium cold stellar matter; b) Y v = 0.5, c) Y p — 0.3. The thick lines refer to the Thomas-Fermi (TF) calculation and the 
thin red lines to the coexisting-phases (CP) calculation. 
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FIG. 2: Transition density, for several temperatures T = 0, 5, 7 and 8 MeV and proton fractions Y p = 0.5 (upper plots), 
0.3 (lower plots) and /3-equilibrium (eq-/3) at T=0 (left part of upper plots), obtained using different methods within the 
models a) NL3 and b) TW. The methods used are: estimation from binodal (Bin), estimation from dynamical spinodal (Sp-d) 
estimation from thermodynamical spinodal (Sp-th), calculation using the coexisting- phases method (CP) and calculation using 
a Thomas- Fermi method (TF). 
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TABLE I: Pasta phases within NL3 and TW for several temperatures and isospin asymmetries. The onset densities of the 
various phases (droplets, rods, slabs, tubes and bubbles) are shown together with the pressure of the homogeneous phase at 
the lower (P 1 ) and upper (P 2 ) border of the pasta phase. The results have been obtained within the coexisting-phases method 
(CP), corrected according to the text explanation, and the Thomas- Fermi approximation (TF), coincident with results in Phys. 
Rev. C 82, 055807 (except for two misprints for P 2 at T = 5 MeV and with Y t = 0.4). 
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